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ABSTRACT 

Tools have been developed to model and simulate the effects of lunar landing vehicles on the 
lunar environment (Metzger, 2011), mostly addressing the effects of regolith erosion by rocket 
plumes and the fate of the ejected lunar soil particles (Metzger, 2010). These tools are being 
applied at KSC to predict ejecta from the upcoming Google Lunar X-Prize Landers and how 
they may damage the historic Apollo landing sites. The emerging interest in lunar mining poses 
a threat of contamination to pristine craters at the lunar poles, which act as “cold traps” for wa- 
ter and may harbor other valuable minerals Crider and Vondrak (2002). The KSC Granular 
Mechanics and Regolith Operations Lab tools have been expanded to address the probability 
for contamination of these pristine “cold trap” craters. 

The KSC tools include simulations of 3D flux of rocket plume ejecta. The trajectories of such 
ejecta can be mapped onto cold trap craters to predict deposition for expected lunar landings, 
assuming that the collection efficiency of the 40K cold trap surfaces is 100%. The processes 
addressed by the KSC tools have now been expanded to address volatile contamination of the 
lunar surface (Stem, 1999). Landing nearby such a crater will result in the migration of signifi- 
cant exhaust plume gas into the cold trap of the crater, and will also create an unnatural atmos- 
phere over the volatile reservoirs that are to be studied. Our calculations address: 1) the time for 
the plume-induced local atmosphere above cold traps to decay to normal levels, 2) the efficien- 
cy of gas migration into a permanently shadowed crater when the landing is outside it but near- 
by, and 3) reduction on contamination afforded by moving the landing site further from the 
crater or by topographically shielding the crater from the direct flux of a lander’s ground jet. 
We also address plume volatiles adsorbed onto and driven inside soil ejecta particles from their 
residence in the high pressure stagnation region of the engine exhaust plume, and how their 
mechanical dispersal across the lunar surface contributes to the induced atmosphere. One addi- 
tional question is whether the collection of soil ejecta along the base of a topographic feature 
will produce a measurable plume volatile release distinct from the background. We mostly ad- 
dress item 2). Item 3) is obvious from our results excepting that the removal distances may be 
large, but changes to landing strategy can improve the situation. 


INTRODUCTION 


The fate of small particles and volatiles (e.g. molecules) emitted from the lunar surface are es- 
timated assuming ballistic trajectories following Keplerian orbits. Under ideal conditions, we 
find that such particles either escape from the lunar exosphere, or collide with the lunar surface 
within one orbit about the lunar center of mass. Following Stem (1999), those ideal small parti- 
cles which collide with the lunar surface are assumed to be reattached to the surface until they 
are “re-emitted” into the lunar exosphere by either a) thermal desorption or b) sputtering by so- 
lar photons or charged particles striking the lunar surface. 

The particle velocity distribution is a function of temperature and particle mass, and we assume 
the emitted particles follow the Maxwell-Boltzmann velocity distribution. Using the Jean’s 
Escape Velocity for lunar surface, the fraction of particles which escape the lunar exosphere is 
dependent on the scale parameter of the particle ensemble (kT/m) 12 or (2RT/M) 1/2 , where T is 
kinetic temperature [K], m is particle mass [kg], M is molar mass [g/mole], and k & R are the 
appropriate ideal gas constants for the chosen units. Given additional information on the ex- 
haust rate of release [kg/sec] and total mass of the emitted particles, we calculate the time de- 
pendent concentration of particles in the lunar exosphere, and characteristic time for deposition 
of the “emitted” particles to the lunar surface. Such “in-flight” concentrations are needed to es- 
timate the rate of particle loss due to exospheric collisions with solar UV radiation or the solar 
wind plasma. These calculations are made for ideal conditions assuming a spherical Moon 
with no surface terrain, a spherical gravity field with all mass concentrated at the Lunar Center, 
and no aspherical/asymmetric gravitational effects. The lunar exosphere is considered to be suf- 
ficiently rarified so that there are no collisions after each particle has been emitted from the lu- 
nar surface. The effects of electric fields are also neglected in this paper, so we are addressing 
only the propagation of neutral (uncharged) particles. 

BALLISTIC TRAJECTORIES IN THE LUNAR EXOSPHERE 

Consider a particle or other small mass released from the lunar surface into the lunar exosphere, 
with no limit to its mean free path. Assuming a spherical moon with all mass concentrated at its 
center of mass and no forces other than gravity, this particle will follow a Keplerian elliptical 
orbit. The orbital equations provide the distance r of such particles from the center of the larger 
mass (not center of ellipse) 

r = p (1 + s cos 0) _1 Eq. (1) 

where p is some constant (defined later), and 

s eccentricity, s = (1 - b 2 /a 2 ) 1/2 

0 True Anomaly 

a semi-major axis of ellipse (defined along apoapsis) 

b semi-minor axis of ellipse. 

Assuming that the particle will reach its highest altitude above the lunar surface at the orbital 
apoapsis along the semi -major axis of the elliptical trajectory, then each orbit can be deter- 
mined from initial conditions v 0 and cp 0 at point p as shown in Figure 1 . 
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Figure 1. Elliptical orbital trajectories for particle emitted from lunar surface at point p, 
reaching highest altitude above lunar surface along the semi-major axis (apoapsis), and 
returning to the lunar surface at point q. 

Given initial conditions for particle speed v 0 and direction cp 0 at the lunar surface, then the nor- 
malized angular momentum h for any particle is conserved for a lossless orbit such that 

h = r 0 ' v 0 ‘ cos cp 0 

where h is constant. Note that r 0 = r moon for initial conditions at lunar surface with r moon being 
the radius of an ideal spherical Moon. Angular momentum h ' m is also conserved since parti- 
cle mass m is constant. Following Vinti (1998), the orbit equation defines normalized total en- 
ergy: 

W = !/ 2 v 2 - p / r <0 

where p = G M, G is the Universal Gravitation Constant and M is Total Lunar Mass. The 
semi-major axis a, orbit constant p and eccentricity s are 

a = - p / 2 W 
p = h 2 / p 

s = ( 1 - p / a) 1/2 = ( 1 - v 2 r 2 cos 2 cp / p a ) 1/2 

where Jean’s Escape Velocity for particles leaving the lunar surface and W = 0 is 

Vesc — ( 2 p / rfyioon) Eq. (2) 

Simplifying by substitution of v esc from Eq. 2, then given v 0 and cp 0 at r 0 

r 0 / a = 2 ( 1 - v 0 2 / Vesc 2 ) = f (v 0 ) 

s = [ 1 - 2 (r 0 / a) • cos 2 cp 0 ' v 0 2 / v esc 2 ] 1/2 = f (v 0 , cp 0 ) 

0 = arccos [s' 1 (2cos 2 cp 0 ' v 0 2 / Vesc 2 - 1) ] =f(v 0 , cp 0 ) Eq- (3) 

The length between focus gi and geometric center C in Figure 1 is as, and the height z a of the 
particle above the lunar surface at apoapsis is 

z a = a(l+s)-r 0 > 0 


Particles emitted from the lunar surface with no additional collisions will escape the lunar envi- 
ronment if their velocity exceeds v esc = 2.375 km s' 1 . Substituting u = v 0 /v eS c in Eq. (3), then 

p = 2 u 2 cos 2 (p 0 

where 0 < 0 < rc/2 for p < 1, 0 = n/2 for p = 1, and rc/2 < 0 < n for p > 1. True Anomaly is 
shown as a function of u and cp 0 in Figure 2. This figure shows which fraction of the particles 
are lost to gravitational escape (u > 1), which land in the hemisphere of point p (p < 1) and 
which land in the opposing hemisphere (p > 1). Orbit eccentricity is shown in Figure 3. The 
trivial case where cp 0 = 0 and u 0 = 2~' A defines a circular orbit following an ideal path along the 
surface of a spherical Moon. Although theoretically possible, such trajectories are not likely 
due to terrain interception, and orbital variations due to asymmetric gravitational fields. 


We developed an Excel spreadsheet and FORTRAN algorithms to calculate semi-major axis a 
[km], eccentricity s, height above lunar surface at apoapsis z a [km], Great Circle path 20 [radi- 
ans], and time in flight t [sec] between initial emission point p and final deposition point q. 
Sample results are shown in Figure 4. The heights of the orbital apoapsis Above Lunar Surface 
(ALS) are shown in Figure 5. 



■10 20 

■ » <0 
□44- M 

□tt-M 

□»-« 

■ *> ■ un 
■IMMI4 
0114- 1# 

■I1B-34D 

■ 140-1U 

■ l»]& 

■ J60 1TO 
DlW-lB 


Normalized Velocity u* = V D / v m „ 


Figure 2. The velocity-elevation partition 
of True Anomaly 0 for particles emitted 
from point p on the lunar surface. 



Figure 4. Particle trajectories assuming 
ballistic elliptic orbits for <p 0 = 30° and 
selected initial velocities. 



Figure 3. Ballistic orbit eccentricity as a 
function of normalized velocity u 0 and azi- 
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Figure 5. Apoapsis height in meters Above 
Lunar Surface (ALS) as a function of nor- 
malized velocity u 0 and azimuth angle <p 0 . 


THE MAXWELL-BOLTZMANN VELOCITY DISTRIBUTION 


Assuming a Maxwell distribution, then the statistical ensemble of particle velocities will be a 
function of kinetic temperature T [K] and particle mass m [kg], such that: 

B(v) = 4/Vti ' v 2 / v p 2 ‘ exp (- v 2 / v p 2 ) 

2 2 1 

where k is the Boltzmann constant [kg m s' K' ]. The peak value or “characteristic scale pa- 
rameter” for this distribution is 

v p = ( 2 k T / m ) 1/2 = ( 2 R T / M ) 1/2 

for ideal gas constant R = 8314 m 2 s' 1 K' 1 and molecular mass for the ideal gas M [g mole' 1 ]. 
The list of mostly neutral particles is derived from the expected products of a popular hypergol- 
ic fuel, which can be found in Table 1. The fraction of particles exceeding the lunar escape ve- 
locity v esc for characteristic molecules and temperatures is given in Table 2. 


Table 1 - Characteristic Exhaust Products of a Common Hypergolic Engine Propellant 


Propellant A 

Propellant B 
(Oxidizer) 

Exhaust Products (Haas, 1984) 

Aerozine 50 
(50% UDMH, 
(CH 3 ) 2 N(NH 2 ) 50% 
Hydrazine, N2H4) 

Dinitrogen Tetroxide 
(NTO), N204 

H20, N2, CO, C02, H2, OH 
Note: relative concentrations of combustion 
products appears to depend on fuel reaction 
temperatures. 


Given molar mass M and kinetic temperature T, the scale velocity parameter v p is used to cal- 
culate the velocity probability distribution B(v). The cumulative probability of B(v) is used to 
estimate the fraction of the distribution exceeding a characteristic velocity such as the Jean’s 
Escape value for surface emissions. Sample B(v) distributions are shown in Figures 6 and 7. 


0.004 

2H1 016 

0,0035 



Figure 6. Maxwell Boltzmann distribution 
for 2H1 016 (18 g mole" 1 ) velocity at vari- 
ous temperatures, including daytime lunar 
surface (380K). The lunar escape velocity 
is shown as the dashed line to the right. 
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Figure 7. Maxwell Boltzmann distribution 
for velocities of Table 1 exhaust products 
at daytime lunar surface temperature 
(380K). The lunar escape velocity is shown 
as the dashed line to the right. 


Table 2. Probability of escape for selected volatiles in Lunar Exosphere 


particle 

M [g mole' 1 ] 

TfKl 

v P [m s' 1 ! 

Prob (v > Vesc) 

2H1 (H 2 ) 

2 

120 K 

998.8 

9.62E-01 % 

380 K 

1,777 

30.7 % 

600 K 

2,234 

51.5 % 

1200 K 

3,159 

76.7 % 

1800 K 

3,868 

85.9 % 

HI 016' (OH) 

17 

120 K 

342.6 

none 

380 K 

609.7 

9.93E-05 % 

600 K 

766.1 

2.23E-02 % 

1200 K 

1,083 

2.12% 

1800 K 

1,327 

9.06 % 

2H1016 

(H 2 0) 

18 

120 K 

333.0 

none 

380 K 

592.5 

4.14E-05 % 

600 K 

744.5 

1 .29E-02 % 

1200 K 

1,052 

1.63% 

1800 K 

1,289 

7.66 % 

2N14or C12 
016 

(N 2 or CO) 

28 

120 K 

267.0 

none 

380 K 

475.0 

6.13E-09 

600 K 

596.9 

5.22E-05 % 

1200 K 

844.2 

1.13E-01 % 

1800 K 

1,034 

1.37% 

C12 2016 
(C0 2 ) 

44 

120 K 

213.0 

none 

380 K 

379.0 

none 

600 K 

476.2 

6.90E-09 % 

1200 K 

673.4 

1.44E-03 % 

1800 K 

824.8 

7.91E-02 % 


TIME OF FLIGHT FOR LUNAR EJECTA 

We derived two time-of-flight formulations for particles leaving the Lunar surface with veloci- 
ties v 0 not exceeding lunar escape velocity v esc , traveling along orbits with constant total ener- 
gy, and returning to the lunar surface without collisions or other interactions along the way. 
The geometry of such orbits is shown in Figures 8 and 9. Kepler’s Second Law is applied, 
wherein the area dA/dt swept out by the vector r over the True Anomaly 0 in equal times is a 
constant. The total area of an ellipse with semi-major axis a and semi-minor axis b is 7iab. This 
total area corresponds to the orbital time period for a particle returning through periapsis to 
point p: 

T = ttT 0 (1-u 0 2 ) - 3/2 

where T 0 = r 0 / v esc ~ 12.16 min, and u 0 = v 0 /v eS c is a normalized particle velocity. Using frac- 
tional area f = 2(Ai + A 2 )/7iab as shown in Figure 8, then the travel time from p to q is: 

T pq =f T 0 (l-u 0 2 ) - 3/2 


; o < f < n 


v a 



Figure 8 - Ideal lunar orbital trajectory 
for particle leaving lunar surface at point 
p, and moving through the apoapsis and 
returning to the lunar surface at point q. 



Figure 9 - Geometry for particle travel- 
ing along a partial orbit from pi to P 2 . 
Whenvectors ri and r 2 are known, then a 
family of orbits will satisfy 0 2 - 0i = const. 


where 

f = arcsin (s) + s (1 - s 2 ) /2 - r 0 2 /ab ’ cos 0 sin 0 + Vi n Eq. (4) 

s = r 0 /a ‘ cos 0 - | s | 


and - 1 < s < 1. Note that f = n when 0 = n and s = - 1 (full orbit from periapsis), and f = 0 
when 0 = 0 and s = 1 (short hop at apoapsis). We also derive the time of flight using the orbit 
equation (Eq. 1) using the geometry shown in Figure 9, which has the following solution 


p 2 rn d9 
ab ^ 0 (l+g cos#) 2 


with 


p_ 

ab 


b 3 /a 3 


Integrating from 0 to n, then: 

f = 2 arctan tan-) — 2s Vl — £ 2 (— — — ) Eq. (5) 

Wl+E 2/ V1+8COS0/ 'l v / 

where 

lim^ „ 2 arctan tan^ ) = n 


The time of flight from any True Anomaly 0 to apoapsis is proportional to f and can be derived 
from either Eq. (4) or Eq. (5). Given ri and r 2 and the Great Circle angular distance 80 i 2 = 0 2 - 
0i between them, then there is a family of elliptical orbits which will pass through both pi and 
p 2 . Each of these orbits has a unique eccentricity s, so the apoapsis height a(l+|s|) and the time 
of flight f(s) for each of these orbits are also unique functions of s. For this analysis, we limit 
the solutions to rl = r2 and confine the source to the lunar surface. 


We solve for f(s) given the case were ri = r 2 = r 0 . In such cases the apoapsis of any orbit will be 
located at 50 i 2 /2, such that 

a = r 0 (1 - s 2 )' 1 { 1 - 1 s | cos (80 i 2 12) ; 0 < | e | < 1 Eq. (6) 

We now have the parameters needed to calculate f(e) for each s, namely 0 = 80 i 2 /2 where 0 < | 
s | < 1, a from Eq. (6), and f from Eq. (4) or Eq. (5). The time of flight f(s) is therefore a 


monotonic function of 8, with the shortest time defined by s = 0 (circular trajectory from pi to 
P2). The longest time is obtained as 8 approaches 1 . Each orbit also defines a unique cp 0 and v 0 
at pi. The fate of volatiles in the lunar exosphere can now be obtained by convoluting the time 
of flight with the velocity distribution as a function of source temperature. We simplify the 
problem by assuming that emissions from the lunar surface are isotropic (uniform in angle). 

Our FORTRAN program EMISSION.F90 integrates the deposition from a point source to 
Great Circle arc lengths 20, summing up all contributions over velocity from 0 to v esc and all 
emission elevation angles cp from 0 to tt/ 2. Results for H20 gas emissions from a point source 
are shown in Table 3 for selected gas temperatures. Contours for such “first hop” deposition 
from a point source are visualized with Google Moon in Figure 10 for water at the characteris- 
tic lunar daytime surface temperature of 380 K. It is the inverse problem that we need to solve, 
namely for each 0 find the relative contribution of emissions integrated over the velocity distri- 
bution. At this point, such a solution is found numerically, with results shown in Figure 10. 


Table 3. H20 | 

gas deposition from point source as function of Great Circle arc [degrees] 

Source 

Temperature 

Fraction lost to 
space 

50% Deposition 
after one hop 

67% deposition 
after one hop 

98% deposition 
after one hop 

120 K 

none 

O 

^1 

O 

1.2° 

4.6° 

380 K 

4.14E-05 % 

2.3° 

3.9° 

18.7° 

600 K 

1.29E-02 % 

U> 

bo 

0 

6.7° 

39.9° 

1200 K 

1.63 % 

OC 

ON 

0 

16.0° 

139° 

1800 K 

7.66 % 

13.4° 

25.8° 

164° 



Figure 10. Deposition contours from a point source in Google Moon, using the results for 
non-ionized H20 gas at 380 K. Trajectories (white) are shown for one hop deposition of 
H20 volatiles at v p for 380K, u 0 = 0.256. Contours are shown as range rings (red) includ- 
ing 50% and 67% of the deposition, from Table 3. 


GEOMETRY FOR LUNAR DEPOSITION 


The deposition of engine exhaust products is estimated for a lunar lander descending to the sur- 
face along a ground track s with engine height h and engine axis tilted down from the local hor- 
izontal by angle j3 0 - The density of engine exhaust is provided by CFD case studies for an 
Apollo Lunar Lander (LM), and these results are combined with an ideal spherical geometry 
model for intersection of rays from the engine exhaust port onto the surface at point p defined 
by along track Si and cross track n t coordinates. 

CFD Case Study - Using CFD simulation for Apollo LM descent engine operating at 45m 
above the Lunar surface, results have been interpolated to a regular grid by 3 -point linear trian- 
gulation. Results are shown in Figure 1 1 for the LM descent engine parameters, with an arbi- 
trary color scheme to highlight the general structure of the derived fields (each color scheme is 
unique). The combination of density and radial velocity at 25m ALS (Above Lunar Surface) 
provides an estimate of engine exhaust flux in units of kg m 2 sec’ 1 which is shown in Figure 13. 
This is a compromise, since the 25m ALS height is the lowest altitude in this simulation which 
provides a cross section of the entire exhaust plume, but may be far enough from the engine to 
provide a “far field” estimate of that flux. We note that deposition calculations will use dis- 
tances from the engine on the order of 10 km to 100 km. 


Figure 12. Engine Exhaust flux [kg m 2 
sec" 1 ] as a function of off-axis angle 0 at a 
distance of 20 m. Approximately 90% of 
the exhaust is contained within off-axis 
angle of 21.2°, and 66% within 12.7°. 

Deposition Geometry - Consider an engine E located at height h above point s 0 , tilted down- 
ward by angle p o so that the exhaust axis is aligned along ray E E’ and oriented above the 
along-track coordinate s, as shown in Figure 13. Following are equations which estimate the 
off-axis angle 0 for point p on the Lunar surface, using the natural coordinates {si, ni} for p 
where s is along-track distance [m] and n is cross-track distance [m]. The following derivation 
is designed to calculate deposition to a rectangular grid in natural coordinates {s, n} for an ex- 
haust source located at height h above grid point {s 0 , 0}. Using each point p defined by natu- 
ral coordinates {si, ni}, the the Great Circle arc v|/ from s 0 to p is 

cos \| / = cos (sfR) cos (ni/R) ; v|/ > 0 

where R is Lunar radius. Note that \\ / defines a range ring about point s 0 with radius s = \|/R, 
which is the dashed line ps’ in Figure 13. Point p is at range r from exhaust source E, where 



pressure density axial v*l radial vel temperature 



engine at h = 45m and descending vertically 
to lunar surface. Density and Axial Velocity 
at height 25 m ALS are used to estimate the 
far-field engine exhaust flux. 


r/R = [1 + (1+h/R) 2 -2(1+ h/R)cos \|/ ] 1/2 


The tilt of this conical surface downward from horizontal is defined by angle p, where 

cos p = (r/R)' 1 sin \|/ ; 0 < p < nil 

Note that r defines a cone with apex at source E and intersecting the Lunar surface along arc 
ps ’ as r is rotated about the z axis through angle y, where 

sin y = sin (ni/R) / sin \\i ; 0 < y < 71 

The off-axis angle for point p is then 

cos 0 = sin p sin p o + cos p cos p o cos y 

The exhaust flux E(0) to point p is derived using the approximation 

SE(0) = 16 kg sr' 1 sec' 1 exp ( - 8.886 rad' 2 0 2 ) 8Q 

and the solid angle subtended on the surface is 

5Q = 5s 8n sin p / r 2 

There is a limitation for rays intersecting the lunar surface, which requires that angle a for any 
point p along the arc ps’ must be greater than nil. Otherwise the exhaust particles will miss 
the surface and either go into orbit or escape the lunar environment. This limitation is defined 
by 

a = nil + p - \|/ ; a > nil or p > \|/ 



Figure 13. Spherical geometry for 
engine exhaust axis tilted down- 
ward by angle p o and axially sym- 
metric about ray E E’ . The off-axis 
angle 0 is estimated for point p at 
natural coordinates {si, ni}. 


Deposition for Apollo 17 Lunar Module Descent - The deposition rate [km m-2 sec-1] to the 
lunar surface is calculated using the known descent pattern and engine parameters for the Apol- 
lo 17 Lunar Module (LM). The instantaneous deposition pattern to the surface at start of brak- 
ing is shown in Figure 14. The LM will land at 720 sec and 522 km down range after start of 
braking. Using the LM descent profile shown in Table 4, the instantaneous deposition pattern 
assuming bursts every 120 sec is calculated and shown in Figure 15. These results are shown 
combined with the LM descent geometry and exhaust cones using Google Moon in Figure 16. 
The full deposition map is shown on normal coordinates in Figure 17, and is shown remapped 
to Google Moon in Figure 18. 


Table 4. Apollo 17 LM Descent Parameters (used in Figure 16) 


time 

s [km] range from braking 

beta [deg] engine tilt 

h [km] engine height 

0 sec 

0 

0 

17.2 

120 sec 

190 

12.0 

14.8 

240 sec 

401 

12.6 

12.6 

360 sec 

446 

22.0 

9.5 

480 sec 

502 

28.9 

6.1 

600 sec 

520 

59.0 

0.86 

720 sec 

522 

90.0 

0 


E 


180 



360 


Along Track Range [km] 


I I below threshold 

■ 0.00001 - 0.201872778 

■ 0.201872778 - 0.302809167 

■ 0.302809167 - 0.403745556 

□ 0.403745556 - 0.504681945 
O 0.504681945 - 0.605618334 

□ 0.605618334 - 0.706554723 
.Q. 0 : 7 .°®554723 _0.807491112 
"□ 0507491112 - 0 .908427501 

□ 0.908427501 - 1.00936389 

□ 1.009363891 - 1.110300279 

□ 1.11030028 - 1.211236668 

□ 1.211236669 - 1.312173057 

□ 1.312173058 - 1.413109446 

■ 1.413109447 - 1.514045835 

■ 1.514045836 - 1.614982224 

■ 1.614982225 - 1.715918612 

■ 1.715918613 - 1.816855001 

■ 1.816855002 - 1.91779139 

■ 1.917791391 - 2.018727779 


0) 

HO 

c 


180 



Figure 14. Apollo 17 LM de- 
scent engine deposition rate 
[g km' 2 sec' 1 ] for 1 sec at start 
of braking at t= 0 sec (LM at 
upper left). The distance cross 
track is shown from 0 (top) to 
180 km (bottom). 


Figure 15. Apollo 17 LM de- 
scent engine deposition rate 
[g km' 2 sec' 1 ] at selected times 
(every 120 sec) with braking 
starting at t= 0 sec (upper 
left). Range shown along track 
is 720 km, with braking start- 
ing at s = 0 km (left) and land- 
ing at s = 522 km. 



Figure 16. Lunar deposition 
bursts from Figure 15 re- 
mapped onto Google Moon, 
with LM positions and model 
of exhaust jet cone. 
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Figure 17. Apollo 17 LM descent engine total deposition [g km' 2 ] after braking starting at 
s = 0 km (left) and landing at s = 522 km. Range shown along track is 720 km, The dis- 
tance cross track is shown from 0 (left) to 180 km (right). 



Figure 18. Total Lunar deposition from Figure 17 remapped onto Google Moon, without 
engine exhaust cones. 



SUMMARY AND RECOMMENDATION 


The techniques demonstrated in this paper show how deposition maps can be constructed from 
engine descent plans. These deposition maps make it straightforward to estimate direct exhaust 
deposition to cold traps. Additional techniques support the calculation of secondary deposition 
after “hopping”, provided that the thermal characteristic temperature at each hop is known. The 
temperature and number of hops enables an estimate of loss during hopping, which is mostly 
due to photo-ionization, which has a time constant at 1 AU ~1.2 day (Stem, 1999). We are cur- 
rently researching the characteristic time constant for thermal desorption, but this time constant 
is likely less than 1 hr which may not be significant to estimates of migration to cold traps. 

These results inform potential future missions where contamination by exhaust gases is not de- 
sired. We recommend that most braking be performed while the braking engine is pointed over 
the horizon, with exhaust gas velocity » lunar escape velocity. This will ensure that most ex- 
haust leaves the lunar environment. Once the hypothetical lander has lost sufficient momentum, 
the braking engine apparatus and excess fuel could be dumped to reduce landing mass. This 
assumes that the braking engine will not be reused for later ascent. It is also important during 
any masss dumping operation that any remaining excess fuel does not escape. The lander can 
then descend with a lower mass (less the descent engine and excess fuel) and minimize exhaust 
during the final descent to the lunar surface. 
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